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Abstract 

The flavor sector of Little Higgs models based on product groups, notably the Littlest 
Higgs with T parity (LHT), has been extensively studied and some amount of fine tuning 
was found to be required to meet the experimental constraints. However, no such attention 
has been paid to other classes of models. Here we analyze the phenomenology of flavor 
mixing in the lepton sector of a simple group model, the Simplest Little Higgs (SLH) . We 
obtain the Feynman rules of the SLH in the 't Hooft-Feynman gauge up to the necessary 
order and calculate the leading contributions to the rare processes fi — > ej, fi — eee and 
fi — e conversion in nuclei. We find results comparable to those of the LHT model, because 
in both cases they arise at the one-loop level. These require the fiavor alignment of the 
Yukawa couplings of light and heavy leptons at the per cent level or an effective scale of 
around 10 TeV. 
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A Loop functions 



1 Introduction 



It is well known that the Standard Model (SM) Higgs mass receives quadratically divergent 
radiative corrections dependent on the cutoff scale of the model. Since electroweak precision 
data (EWPD) requires that the Higgs mass be at most of the order of the electroweak scale 
V ~ 246 GeV, naturalness arguments demand that this cutoff scale be near the TeV to avoid 
fine tuning the model parameters (hierarchy problem). Hence, new physics effects are expected 
before or at the TeV scale. Many beyond SM scenarios have been proposed to solve this 
hierarchy problem, such as super symmetry, technicolor and extra-dimensions among others. 
However, the same EWPD and flavor physics in general disfavor new particles at scales 
somewhat below ~ 10 TeV giving rise to the so called little (flavor) hierarchy problem. 

Little Higgs (LH) models [IHl] are also an attempt to solve the hierarchy problem, bridging 
the gap to ~ 10 TeV. This is done making the Higgs a pseudo-Goldstone boson of a new 
approximate global symmetry broken at a scale / ~ 1 TeV. The Higgs mass is protected from 
the one-loop quadratically divergent corrections and leaves only a two-loop sensitivity to a 
10 TeV cutoff which is not considered unnatural. Beyond this scale an unspecified ultraviolet 
completion is needed but this completion can now elude EWPD and flavor constraints. 

There are several possible realizations of LH models depending on how the new symmetries 
are implemented. Broadly speaking, they can be separated into two categories: those that 
have the SM gauge group emerging from the diagonal breaking of the product of several groups 
(for instance {SU{2) x U (1))'^) and those where it emerges from the breaking of a larger simple 
group (for example SU{N) x U{1)). The former are called product group models while the 
latter are termed simple group models. There are some general features characteristic of each 
of the approaches [5]. In the absence of further symmetries, product group models have 
additional free parameters from the gauge couplings of the different groups that diagonally 
break down to the SM gauge group, the SM Higgs can be embedded in definite models into a 
single sigma-model multiplet and it is possible to make the fermion sector relatively simple. 
On the other hand, simple group models have the enlarged gauge group couplings fixed by the 
SM values but they require, in specific models, at least two sigma-model multiplets. Moreover, 
fermion multiplets must be extended to transform under the new gauge group. 

At any rate, LH models introduce new particles with masses of order /. However, EWPD 
generally requires / > 4 TeV [SHS], reintroducing the little hierarchy problem. In product 
group models, these constraints can be alleviated by introducing an additional discrete sym- 
metry, T-parity [10H12| . This symmetry exchanges the gauge groups, making almost all new 
particles T-odd and all SM particles T-even. In this way, dangerous tree-level couplings of 
light fields with only one heavy particle are forbidden. This avoids large contributions from 
higher dimensional operators obtained by integrating out the new heavy fields, thus relaxing 
the tension between EWPD and a lower /, which can be of the order of the TeV in this 
case [13]. As a byproduct, the predicted (now relatively light) new particles could be eventu- 
ally observed at high-energy colliders [Ulll]. This justifies the attention which these models 
have received. In contrast, simple group models have no consistent way of introducing a 
similar mechanism [12j, and thus they are somewhat disfavored. 

However, one must also consider fiavor constraints which are in general more stringent 
and translate into more restrictive, and also complementary, limits on the model parameters 
(see [IS] for a review and further references). Flavor violating processes depend on the new 
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heavy scale as well as on the misalignment of SM and heavy flavors. Hence, the corresponding 
limits can be satisfied sending / to a high enough value or aligning both sectors with a high 
enough precision. For instance, in the Littlest Higgs model with T-parity (LHT) present 
bounds on lepton flavor violating (LFV) processes require / > 10 TeV or a misalignment 
of at most 1% between the SM and the heavy fermion mass matrices |16H19j . These limits 
are not stronger because T-parity forbids tree- level flavor changing Z couplings. At any 
rate, flavor changing neutral currents (FCNC) involving only SM external fields are induced 
by the exchange of (heavy) T-odd particles at one loop, resulting in the previous bounds 
and reintroducing a little flavor hierarchy problem. Flavor violation in the quark sector has 
also been addressed in the literature, both in the Littlest Higgs with |20fl22] and without T 
parity [23]. 

Once / is of the order of several TeV, it is of the order of the scale implied by the EWPD 
bounds on simple group models. This means that simple group models and the LHT would 
be on similar footing as long as flavor constraints on the former models are not more stringent 
than in the LHT case. As we shall show in the following, the LFV limits on the Simplest 
Little Higgs (SLH) model j24U25] are comparable to those on the LHT. This is so since in this 
simple group case the matter content of the model guarantees the absence of tree- level charged 
lepton FCNC, and the corresponding LFV processes are then one- loop suppressed. In the 
SLH model the new global symmetry is (S'C/(3) x where only the diagonal subgroup 

SU{3) X U{1) is gauged. This gauge symmetry is broken at the scale / into the SM gauge 
group SU{2)l X U{1)y- Left-handed (right-handed) matter flelds transform as SU{3) triplets 
(singlets) , implying only the addition of heavy quasi-Dirac neutrinos to complete the lepton 
multiplets |261l27j. Hence, only neutrinos have tree-level FCNC, although mixing only light 
and heavy neutrinos at the order considered, with no immediate observable eff'ect [28l[2?j . 

The quark sector in the SLH is more involved. There are two ways of embedding the 
SM quark doublets into the new SU(3) triplets. In any case, the quark Yukawa Lagrangians 
allow for mixing between heavy and light quarks of all three families. Since in this paper 
we concentrate on the contributions to the basic LFV processes — )• e7, — )• eee and //~e 
conversion in nuclei, the quark sector is only relevant to the last process. For simplicity, we 
will suppress all quark related mixings in this case to evaluate only the lepton mixing effects. 

We flnd that predictions for the SLH model are similar to those of the LHT. For instance, 
the present bound on Au — t- e Au requires / > 14 (16) TeV for the anomaly-free (universal) 
quark embedding and natural values of the other model parameters in the SLH case, to be 
compared with 10 TeV for the LHT model. Alternatively, the misalignment parameter must 
be sin 20 < 0.005 (0.004) in the former case, and sin 20 < 0.01 for the latter. The processes 
— > e7 and fi — )■ eee give comparable but less stringent limits. LFV r decays will not 
be presented here, as they are less restrictive, since the branching ratios are suppressed by 
a factor B{t — )• ii'ei'r) ~ 0.2 and the experimental limits are several orders of magnitude 
weaker. 

LH models, as technicolor models in the past [30], exhibit a naturalness flavor problem, 
more demanding in principle than the little hierarchy problem. Obviously, one can argue 
that they have signiflcative regions of parameter space allowed by present flavor constraints, 
but these are characterised by a large scale or a small mixing angle, or a combination of 
them, when using a relevant physical parameterization. Also some observables can have 
cancelling contributions in other regions but these are relatively narrow and not common to 
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all of them. They are prime candidates for explaining any future observation of (lepton) flavor 
violation beyond the SM. Thus, it is important not only to investigate the requirements for 
constructing realistic LH models, but to eventually interpret any possible departure from the 
SM predictions which may be observed at ongoing or planned high precision flavor experiments 
like MEG [SUES] at PSI and PRISM/PRIME [MlIM] at J-PARC. A general discussion of 
different scenarios, including the LHT model as well as supersymmetric and extra dimensional 
models, can be found in ^35u36j . The comparison of muon and r branching fractions, when 
the latter are known with a better precision, should eventually help to discriminate between 
specific models [37j. If one just relies on processes involving only the first two families, models 
can be essentially classified depending on the order where fj, — )• eee appears, i.e. tree level, one 
loop, or beyond, compared to /i — >■ e^y, which is always loop suppressed as required by gauge 
invariance. On the other hand, the use of the polarization of the initial lepton provides further 
observables to determine the model structure, as recently discussed in the LHT case [38]. 

Here we will present limits on the SLH model derived only from LFV effects, which are 
similar to those of the LHT, as already emphasized, because in both cases they are one-loop 
suppressed. In Section [2] we detail the calculation of the necessary Feynman rules in the 
't Hooft-Feynman gauge to obtain the prediction of the LFV processes in the SLH for the 
first time. In Section [3] we present the general structure of the LFV branching ratios and 
form factors and calculate the different contributions of the SLH. Section 2] is dedicated to 
describing the numerical predictions of the model, while our conclusions are summarized in 
Section O 

2 Feynman Rules for the SLH model 

Our objective is to study lepton fiavor violation in the SLH model. The only new source of 
LFV in the SLH model is the misalignment of the SM down-type lepton mass matrix with 
the new heavy neutrino mass matrix. Since we also neglect SM neutrino masses and mixing 
effects, this is in fact the only source of LFV. This means that mixing matrices only appear in 
vertices that couple SM leptons to the new heavy neutrinos and, since our external particles 
are charged (e, fi and quarks), only charged currents can contribute to the flavor change. An 
extension of the SLH model including the observed neutrino masses can be found in [26]. 

The SLH model is introduced in |24p25j and the interaction Lagrangians are derived in [5] . 
However, some additional manipulation is required since we prefer to work in the 't Hooft- 
Feynman gauge. This forces us to keep the would-be-Goldstone bosons (WBGB) explicit when 
obtaining the Feynman rules. Also, the LFV matrix elements are of order v^//^ in general so 
we have the additional complication of obtaining all the couplings up to the necessary order 
to guarantee the consistency of the calculations to this degree of precision. 

Our first step is to obtain all the physically relevant fields in the model. When the fields are 
originally defined, they are the interaction states and usually are either not mass diagonal or 
have non canonical kinetic terms or both. This is especially true when we take our expressions 
to order v'^/f'^. Many of the fields would need no transformations and would directly be the 
physical fields if there were no electroweak symmetry breaking (EWSB). In many cases the 
notation is inspired by this fact. Also, as expected in a simple group model, several constants 
that in principle have no specific value end up being determined by the requirement that the 
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model must contain the SM within it. 



In the following we define our basic interaction fields and analyze in turn the different 
sections of the Lagrangian. 



2.1 Basic fields and expansions 



The SLH model is an SU{3)l x U{1)x gauge theory. This gauge symmetry is a diagonal 
subgroup of a global {SU{3) x U{1))-^ x {SU{3) x U{1))^ group. The global symmetry is 
spontaneously broken to (^SU{2) x U{1)^ ^ x (S'[/(2) x U (1)) ^ and the gauge symmetry reduces 
to the SM gauge group SU{2)l x C/(l)y. We will write the covariant derivative for this gauge 
group in the following form: 



gtw 



(2.1) 



with tw = sw/cw the tangent of the Weinberg angle 9w- Writing the SU{3) generators 
in the fundamental representation (3), the SU{3) part works out as follows: 
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Here we have introduced the definitions of the gauge bosons 

y'^^. Some of these will have to be rewritten in terms of the massive physical gauge fields 
later on. The value of the gauge coupling g^ ends up being set by the requirement that the 
photon couple to the electric charge. 

The scalar sector of the SLH is a non-linear sigma model. There are two scalar multiplets 

$1^2 transforming as (3, 1) and (1,3) under S\J{^^\ x S'?7(3)2, respectively (each then a (3) 
representation of the S\J (3) gauge group with a U (l)x hypercharge of —1/3) that include the 
SM Higgs doublets as well as new Goldstone bosons. They can be expressed as follows: 
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with tp = spjcfi the ratio of the vacuum expectation values of the two Higgs triplets. To 
introduce EWSB we will substitute /i° = + if)/\/2 — ix and = — (/>^. Note that the 
term exp (^^^ can be rotated away using a gauge transformation SU (3)l x [/ (l)x as would be 
the case if we wished to work in the unitary gauge. The structure of the Q and 0' matrices is 
determined by the broken generators of the gauge symmetry. One can add identity matrices 
from the U{1) group which allows us to set the first two elements of the diagonal matrix 
independently from the third element (/t and k,'). Certain choices are convenient [5j but 
they are not needed in our calculation since they only affect neutral Goldstone bosons. The 
scalar triplets are then expanded up to any given order in v/f (when we include EWSB) and 
expansions of up to order w^//^ are required in some cases to obtain the necessary precision 
for the couplings. 

To build the fermion sector of the model, the SM fermions have to be included in repre- 
sentations of the larger SLH gauge group. The simplest way to do this is to embed the SM 
fermions into SU{3) triplets. The fermion sector of the SLH model can then be broken down 
as follows. 



Each lepton family consists of an SU (3) left-handed triplet (3) and 2 right-handed 
singlets (1). There is no right-handed light neutrino: 



1^^ = {l^L,iL,iNL)m, iRm, Nr^ . (2.7) 



• The structure of the quark fields depends on the embedding we select. In the so called 
universal embedding, the quark sector is analogous to the lepton sector. The only 
difference is that we have three right-handed singlets: 

Qm = i'^L,dL,iUL)m, URm , dRm , Ujim ■ (2.8) 

However, this universal fermion sector leads to SU{3) and U{l)x gauge anomalies al- 
though these do not affect the SM gauge group. Since the SLH model is an effective 
theory this is not necessarily a problem because one can add additional fermions at the 
cutoff scale of the theory to cancel the anomaly. On the other hand, one can construct a 
quark sector that is directly anomaly- free with no additional degrees of freedom [39U40] . 
This is known as the anomaly-free embedding and it requires that the first two families 
contain SU (3) left-handed conjugate triplet representations (3) and three right-handed 
singlets. The third family is analogous to the lepton sector: 

Qi = {dL,-UL,iDL) , dR, UR, Dr, (2.9) 
Q2 = {sL,-CL,iSL) , SR, CR, Sr, (2.10) 
Ql = {tL,bL,iTL) , tR, bR, Tr. (2.11) 



The gauge representations and hypercharges for the fermion sector for the different embed- 
dings are summarized in table [H 
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Universal embedding (U) 


Fermion 


Ql,2 




URm, Uum 


dRm 




NRm 


eRm 


Qx charge 


1/3 


1/3 


2/3 


-1/3 


-1/3 





-1 


SU{3) rep. 


3 


3 


1 


1 


3 


1 


1 


Anomaly-free embedding (AF) 


Fermion 


Ql,2 


Qs 




C?i?m, DRm, SRm 


Lm 


NRm 


eRm 


Qx charge 





1/3 


2/3 


-1/3 


-1/3 





-1 


SU{3) rep. 


3 


3 


1 


1 


3 


1 


1 



Table 1: Quantum numbers of the fermion fields for the universal and anomaly-free embed- 
dings. 



2.2 Gauge and Goldstone boson sector 

In this section we analyze the pieces of the Lagrangian that involve only gauge and Goldstone 
bosons. 



2.2.1 Scalar Lagrangian 



From the gauge invariant Lagrangian 



(2.12) 



we can readily obtain the charged gauge boson mass terms. As a first approximation we keep 
terms up to order / p. To this order the charged boson sector is diagonal, 
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(2.14) 

(2.15) 



This level of precision is sufficient everywhere except when obtaining the correct 0{v^/f^) 
couplings for Goldstone bosons. We need to go up to order v'^/f'^ to obtain the corrections to 
the W mass and higher order corrections to the gauge boson eigenstates for these couplings. 
Taking our expansion up to order four, the mass terms then read: 
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(2.16) 
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To obtain the correct physical states we must rotate the original fields as follows: 



x^± 




(2.17) 



Note that the physical states W and X differ from the interaction states only by a term 
of order / f'^ . This difference is irrelevant almost everywhere in our calculation, but is 
important in determining the would-be-Goldstone boson states. The interaction fields W and 
X will be considered equal to the physical fields elsewhere. 



Mw = 
Mx = 



The masses of the physical fields work out as: 
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(2.18) 
(2.19) 



The 0{v^/f^) correction to Mx is neglected since it is unimportant for our calculation. 

The neutral sector is already non-diagonal at order C(u^//^) and requires some more 
work: 
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Diagonalizing this matrix, the masses at order 0{v^ / p) are: 
Ma 
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(2.22) 
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(2.26) 



The first order mixing matrix for gauge bosons is then: 

cw 
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V3 



-sw 

Sw_ 



f2 



w I 
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(2.27) 
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Additionally, the physical Z and Z' states also require the replacements: 

Z' ^ Z'^bzZ, 
Z Z-6zZ', 

where 

{l-t 



Sz 



W) 



t 



w 



Sew 



p 



(2.28) 



(2.29) 



We now need to find the actual Goldstone eigenstates. The ones that appear in the original 
expansion have non-diagonal kinetic terms to order j p and there is mixing of these states 
with the gauge bosons through terms V^d^(\). We need only the charged sector since the 
neutral Goldstone bosons do not contribute to lepton flavor mixing processes. 

The kinetic terms for the charged Goldstone bosons and the Goldstone-gauge mixing terms 
read (interaction fields): 






(2.30) 



(2.31) 



We can combine all the necessary Goldstone transformations into a single pair of equations 
that express the original interaction eigenstates in terms of the final Goldstone states (to order 
f^//^) associated to the W and X bosons and have canonically normalized kinetic terms: 



X 



i± 



^ Ti 1 + 




(2.32) 



The calculation of these states required the use of relation ()2.17p to obtain the j p correc- 
tions. 



Taking both (j2.32p and the gauge boson rotations (j2.27p and (j2.28p into account, the 
relevant Feynman rules can now be obtained. The results to order 0{v'^ / p) are given in 
table m 



2.2.2 Vector Boson Lagrangian 

From the Lagrangian 



z 4 g 



(2.33) 
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±\Mx 


JX^X^ 
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Zx^x^ 






2SWCW 1 '-2swV3-t?y 




±\Mwtw ± '^^zMw ^~f^. 


Zcp^cp'f 




x^X^Z' 


=Fi Z^-^, — ^ ± i5z Mx " 


Z'x^x^ 


1 — _ 

T ^ T r)" 




^2.v^V3-t?v- 


(p^W^Z' 


±iMH/ T iSzMwtw 




-r -pA^V ""W 



Table 2: Vertices [SV^V^] = ieRgt"" and [V^S(pi)5(p2)] = ieG(pi - ps)''- 



vvv 
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VVV 


J 


-yX+X- 


-1 




-1 


zx+x- 


'^w~^w &z /q j-2 
2svi'Ci^ 25^ V ^ 


zw+w- 


sw 
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1 /o +2 x_'^w~''iv 
2svK V ^ ^ 
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Table 3: Vertices [V^(pi)V,(p2)Vp(j>3)] = i^b^"(P2 - PiY + 9'''{P3 - P2Y + g^'ipi - PsT]- 

and using (j2.27p and (j2.28p we obtain, in a straightforward way, the [VVV] couplings we need 
for our study of LFV processes. They are given in tabled 



2.3 Lepton sector 

Here we will treat the pieces involving leptons which include their Yukawa Lagrangian as well 
as their gauge interactions. We will obtain the lepton masses and physical states and their 
couplings to gauge and Goldstone bosons. 



2.3.1 Lepton Yukavi^a sector 

Lepton masses follow from the Yukawa Lagrangian: 



£y 3 iX'^NRrn'^lL^ + '-^iR^eijk^mi^ + h.c. , 



iA^ 



(2.34) 



where the quartic term preserves the global symmetry (L^ transforms as (1,3) under SU (3)i x 
S'?7(3)2) and Xn can be taken diagonal after a proper field redefinition. Firstly we need to 
determine the actual physical states of the leptons. Keeping only the mass terms to order 
C(^^//^) we have the following Lagrangian: 



Cy ^ -f 8/3X1^ 



+ C/3^ArWLn + li.c., (2.35) 
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where 



V2ftp 



C/3 



12/2 











C/3 





(2.36) 



The matrices \m and are not necessarily ahgned. Thus, in the basis where the former 
is diagonal, the latter mixes different light lepton flavors. Denoting the eigenvalues of \i as 
y^., the light lepton masses are given by 



(2.37) 



whereas the left-handed components of the light physical fields are obtained by the replace- 
ment: 



(2.38) 





Furthermore, according to (j2.35p each heavy neutrino is mixed just with the light neutrino of 
the same family. To separate them we rotate only the left-handed sector. To order Oiy"^ / f"^)^ 
the physical states for the neutrinos are given by: 



(2.39) 



Notice that the mixing angle 6y between light doublet and heavy singlet neutrinos is experi- 
mentally constrained to be small |29y 41j. Although the bound is flavor dependent (0.05, 0.03 
and 0.09 at 95% C.L. for z^e, t'^t and Vr respectively) we will assume a typical upper limit 
5i, < 0.05 for illustration purposes. 

Since one can safely consider the SM neutrinos as massless, we have chosen to rotate them 
in the same way as the light charged leptons. Finally, the heavy neutrino masses are: 



mN, = fs^X 



N ■ 



(2.40) 



We now need the Goldstone- lepton couplings. These require expansions of up to order 
0{v'^ / f^) in order to get the couplings to order 0{v'^ / f'^). Goldstone rotations from the scalar 
sector in (j2.32p are also needed. We want only terms that couple one charged Goldstone 
boson (x^ and (j)^) with two fermions. Taking these considerations into account and using 
the physical fermion states (j2.38p and (j2.39p we obtain the couplings in table [H Note that 
several coupligs are zero because we have neglected the SM neutrino masses. 



2.3.2 Lepton-Gauge sector 

>Cf = Ipmi^i'm , tpm = {Lm^Rm, NRm} ■ (2.41) 

The expression for the covariant derivative is in equation (12. ip . We note that the only places 
where f^//^ corrections appear are in the definitions of the physical states of the leptons (in 
and in the ZZ' mixing (in 6z)- We must use (j2.38p and (j2.39p and the definitions of the 
physical gauge bosons to obtain the relevant Feynman rules which are listed in table El 
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SFF 


CL 


CR 




1 Tn AT- / Ji'^ \ A A 
1 '"JVi / 1 Ov\ ^iV 


1 h \2\ yij 


V2sw '^x 2 ) 


V2sw ^x \^ '^^1 ^t- 


x-ijNi 




7 cii \ 

V2SIV MxV 2 ; 




r i ruNi i^ij 


X 1 *J T/«J 













X 1 


X ^il^i 


r 1 










i /'i 

^SM. Mw \ 2 ) 


(j)~£ii'i 


i "^^i (i Si\ 






Table 4: Vertices [SFF] = \e{cLPL + crPr) for the lepton sector. 
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^/2sw\ 2 J 







1 ymi 

^ y/2sw ^ 









^IV X_ «IV 
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r 1 ymi 

2SWCW e 







^-V2s^ 







V2SW V 2 j 





y" 























x2 i T^mi 











23M^c^wV3-t'v^ ^ 2.i^c^ 




Z'viVi 


1-2^^ /i x2C^v.(3-t^)\ . 1 





25wc?^V3-t?^ V ^ 1-2^4' ; ^-2.^-0^ 


Z'NiNi 


1 fi ^;f2 3-t^^\ 
^H^V3-t^ I' 2 J 





Z'NmVi 


-<^'^2ii7\/3-*w^r' 






Table 5: Vertices [VFF] = ieigLPi + 9rPr) for the lepton sector. Notice that y° 7^ 
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2.4 Quark sector 



For the //N — t- eN conversion process we also require some of the quark couphngs. The full 
mixing structure of the quark sector is much more complex than that of the lepton sector 
and, in general, all light quarks mix with other heavy and light quarks from every family. 
However, we are only interested in the mixing effects that are a consequence of mixing in 
the lepton sector so we will neglect most mixing effects of the quark sector. Corrections of 
order f^//^ to vertices are only needed for particles involved in triangle diagrams and, since 
quarks only appear in box diagrams, vj f precision is sufficient. This simplifies the calculation 
of our Feynman rules considerably. However, we do have to analyze both the universal and 
anomaly-free embeddings since they produce different results in general. 



2.4.1 Quark Yukav^a sector 

For the anomaly-free embedding, the basic Yukawa Lagrangian reads: 



m 



+iAf 4„Q^<i>i + iAf 4„g^$2 + i^URme^Jk^f^l'Q'^ , (2.42) 

where n = 1,2; i,j,k = 1,2,3 are SU{3) indices; runs over {dji, SR,bji, Dr, Sr) and URm 
runs over {uR,CR,tR,TR)\ it^^g and are linear combinations of tR and Tr; and 
are linear combinations of dR and Dr for n = 1 and of sr and Sr for n = 2: 



2) 



\fcpd^R, + \fspdl, ^^_ ->^fspd\, + Xfcpdl, 

\fc13d\2 + ^fspd^^ ^ -^fsi3d\i2 + ^fcisdR2 

^{Xffcl + {Xffsl ' ^(Af )^4 + {Xffsl 



We require the collective structure with different right-handed quarks entering in the and 
<1>2 quartic Yukawa couplings. By a proper field redefinition, A^ can be taken diagonal in 
general and, for simplicity and to avoid large quark fiavor changing effects, we also assume 
A2 to be diagonal [HIS^ • 

Before the EWSB we obtain the following masses for the heavy quarks: 

mr = f^{X\fcl + {Xifsl , (2.46) 
mn = f^{Xffcl + {Xffsl, (2.47) 
= /^(Af )'c| + {Xf fsl . (2.48) 
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After the EWSB, the quark mass terms work out as follows to leading order: 



^; f 

I J \mn- 

f, n ^ w((Af)MA£)^ ^; AfAf 

-msSRSL - -j= -.SrSl + -7= , . . SRSL 

^^(Af?4TiAff4 ^ ^(Af74TiAf)^ 

+-^|A^Jii^6L + h.c.. (2.49) 

Since we are interested in lepton flavor mixing we will assume no flavor mixing in the 
quark sector, which might otherwise dilute some of the effects we wish to highlight. We 
essentially set all the A™ and A^" that mix different families or heavy and light quarks to 
zero and all others are fixed by the light quark masses. Only the heavy-light mixing within 
each family remains. We neglect terms proportional to / f'^ and rotate the left-handed fields 
to obtain the physical quark states (heavy quark masses get corrections at this order that will 
be neglected as well): 

Tl ^ Tl + SttL , (2.50) 
tL ^ tL~ , (2.51) 
Dl ^ Dl + SddL , (2.52) 
(II ^ dL- SdDL , (2.53) 
Sl ^ Sl + SsSl , (2.54) 
SL ^ SL- SsSl , (2.55) 



where 



S, ^ ^ipM^m , (2.56) 



2 



(2.57) 
(2.58) 



are complex in general. 

Taking all this into account we get the SM quark masses: 

^ / \11 ^ / \22 ^ f \3 /o cnA 

= -^jK , me = --j=jX^ , m, = --^-X, , (2.59) 
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SFF 


CL 


CR 


x~Du 


1 mo 


1 m„ 


V2sw Mx 


V2sw Mx 


x~du 





X* 1 ™d 



SFF 



(j) du 



CL 



Zsw 



CR 



ZSw 



Table 6: Vertices [SFF] = ie{cLPL + crPr) involving first family quarks in the anomaly-free 
embedding. 



mt 



rrid 



V2 



V Af Af 

-^^{Xffcl + iXffsl 

id2 \d2 



Af A^ 



V2 



c} + (Af ) s 



2 2 

/3 



(2.60) 
(2.61) 
(2.62) 



The 5(i^s,t parameters can be expressed in terms of the quark masses: 

1 



vniQ 



"2^/2/2 c^sp 



cl + e 



1 



f 2^,2 



(2.63) 



where the + (— ) sign stands for g = t {d,s) and e = ±1 depending on the corresponding 
values of Ai and A2. We will use this expression to study the predictions of the model because 
it involves parameters with a more straightforward physical meaning. 

With this information and redefining the Goldstone fields as given in (|2.32p , we can obtain 
the relevant quark-Goldstone boson couplings for our processes. Those involving first family 
quarks are given in table [6j Remember that we have removed all quark flavor changing 
vertices so there is no CKM matrix. 

The situation is similar in the universal embedding although the Yukawa Lagrangian is 
different: 



ymn 



(2.64) 



Here m,n = 1,2,3 are generation indices and i,j,k = 1,2,3 are SU{3) indices; dm runs over 
the down quarks {d, s, b) and are linear combinations of the light and heavy up quarks: 

\un„ „,l I \un „ „,2 
Al C/3Ur„ + A2 Si3Uj 



URn 



''Rn 



""^^4 + (AD^s2 



(2.65) 
(2.66) 



Analogously to the anomaly free case, we assume the collective structure. As before. A" can 
be made diagonal by a field redefinition and A2 is also taken diagonal to avoid large quark 
flavor effects. 
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SFF 


CL 


CR 







S* 1 ™di 

«i V2sw Mx 




1 rnu^ 


1 nid^ 


V2sw 


V2sw ^'^x 



SFF 

(f>'^Uidi 



-JJidi 



CL 



i(5„ 



CR 



V2sw 



1 '"dj 

V2sw 



Table 7: Vertices [SFF] = ie{cLPL + crPr) involving first family quarks in the universal 
embedding. 



The mass terms work out as: 



Cy D -f^{\TYcl + {\rrs'pURnULn + ^ 



spcp 



(a: 



2 J 



\ un 
^2 



vf 



UruULt, 



V2 



URnULn + '^J^'^d ^Ri'^Lj + h.C. 



(2.67) 



We have neglected terms proportional to v"^//. 

We will again ignore all generation mixing terms. This means setting Aj' = ^li^ij- The 
only remaining mixing terms involve the light and heavy up quarks of each generation. The 
following rotation of the left-handed fields is required to obtain diagonal mass terms: 

ULn 

where 



ULn — ^UnULn , 



V SpCp[{Xf 



(AD'] 



V2f (Xrfcl + (Ar)2,2 ■ 
The quark masses to order v/ f are: 



V2 
vf 



{xr?cj, + {\T?sj 



(2.68) 
(2.69) 

(2.70) 

(2.71) 
(2.72) 

(2.73) 



V2A ' ■ 

Here (5u„ can also be written analogously to (|2.63p using the global + sign and with e = ±1 
depending on the values of Ai_2- 

The Feynman rules for quark-Goldstone couplings to order v/ f are given in table [71 



2.4.2 Quark-Gauge sector 

For the anomaly-free embedding the Lagrangian is: 

Cf = QmlB^Qm + URmiP^URm + dRmlB^dRm + T/jI^Tr + DriB'^Dr + SriB^Sr . {2.1 A) 
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VFF 



-ydd 
W-Du 
W-du 

Zuu 

Zdd 
X-Du 
X'du 

Z'uu 



9L 

_2 
3_ 

1 
3 



V2. 

Gcwsw 
i 

3— <^ 



"7 



_2 

3_ 

1 

3 



3civ 



Z'dd 



&sw 



3sT4/\/3^ 



Table 8: Vertices [VFF] = ie^giPL+guPji) for the quark sector in the anomaly-free embedding 
entering in our calculation. 



Here we have taken into account that the first two families are in the anti-fundamental 
representation: 



^3/1 



1 



d.-igAlTa + ig.-B-, , 



D 



du + igx 



1 



d^ + ig^-B^ 



For the universal embedding, the Lagrangian is more symmetric: 

C = Qmi£^^Qrn + URml^URm + dRmlB'^dRm + tjR^lEf^URm , 

where 



d^-igA^Ta + igJ^Bl , 



da + i^a: 



- ) K 



d^ + ig^-B^^ . 



(2.75) 
(2.76) 

(2.77) 

(2.78) 

(2.79) 

(2.80) 
(2.81) 
(2.82) 



The Feynman rules for the gauge-fermion couplings in both embeddings are given in 
tables [8] and [9l We only include the vertices and order relevant for our calculations. 
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VFF 


9L 


9R 


^uu 


2 
3 


2 
3 


^dd 


1 

3 


1 
3 


W+Ud 


r* 1 





W ua 


1 

V2sw 


n 
U 


Zuu 






Zdd 


Gcwsw 


sw 


X+Ud 


i 





X+ud 







Z'uu 


3+tw 


3+t^ 




6sw\/3-t2^ 


Z'dd 







Table 9: Vertices [VFF] = ie{gLPL + 9rPr) for the quark sector in the universal embedding 
entering in our calculation. 

3 Lepton Flavor Violating Processes 
3.1 General structure 

Here we summarize the expressions of the final branching ratios and conversion rates in terms 
of the different form factors we will obtain in the following sections. In general we use the 
notation in p^[T9] . 

3.1.1 fj, ^ ej branching ratio 

The partial width for ^ — t- e7 can be expressed in terms of the dipole form factors of the 
lepton flavor changing one-loop vertex as: 

r(/.^e7) = |m3(|F2,|2 + |F^|2). (3.1) 

The branching ratio is then obtained by dividing by the total muon decay width which can 
be approximated by 

3.1.2 — eee branching ratio 

The treatment of the decay width for the SLH is very similar to the one for LHT. The only 
difference here is that we have additional penguin diagrams due to the exchange of Z' gauge 
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bosons. We define the amphtudes and form factors as follows: 

o2 



A^Tpenguin = q2^{Pi) [Q^l^{AfPL + A^Pr) + m^\a^'^ QM^Pl + A^Pr)] u{p) 
Xu{P2hfiV{p3) - (Pl ^ P2) , 



M 



A^Zpenguin = JpKPl) [I'^iFlPL + FrPr)] u{p) u{p2) b^.{ZlPL + Z^Pr)] vips) 

z 

-{Pl ^P2) , 

z 

-[p\ ^ P2) , 

boxes = e^B^ [u{pi)-f^'PLuip)] [uip2h^^PLv{p3)] 

+e^B^ [u{p^)^^Pru{p)] [u{P2h^PRv{P3)] 
+e^B^ {[u{pi)j''Plu{p)] [uip2h,.PRv{p3)] - {pi ^ P2)} 
+e'^B^{[uipi)j''PRuip)] [uip2h^.PLviP3)] - (pi ^ P2)} 

+e^B^{[u{pi)PLuip)] [u{p2)PLvip3)] - {P1^P2)} 

+e^B^{[u{pi)PRu{p)] [u{p2)Prv{p3)] - {pi ^ P2)} 
+e^Bl {[u{j)i)(j^"' Plu{p)] [u{p2)a^,PLv{p3)] - {pi ^ P2)] 
+e^Bf{[u{pi)a^'^PRu{p)\ [u{p2)a ^,Prv{p3)] - {pi ^ P2)} , 



(3.3) 



Fl 



FlIQ 
-Fl, 



2 aR 
> ^1 



Fl/Q' A^ = -{Fl + \Fl)/m, , A^ = -{Fl, - \Fl)/m^ 



Fr 



F'r 



F'r 



'R 



(3.4) 



We can then use the following expression |42p43j to obtain the partial width: 



r(/i 



eee 



where 



9 ^ 
am,, 

327r 



|Af |2 + |ylf 1^ - 2{AiA^* + A^A^* + h.c. 



\Rl2 



L aR* 



\L aR* 



V 3 nie 3 y 



i(|i?f|2 + |Sm + i(|B,^P + |B 



i?,|2\ 
2 I J 



+ ^{\B:^\'+\B^r)+6{\B^\' + \Bi 



Rl2\ 



1 



-{BtBt + B^Bt+h.c.) 



tR T3R* 



+ ^{AfB^* + Pf* + A^B^* + AfB^* + h.c.) 



- -(A^-Bf* + A^B'^* + A^B^* + A^B^* + h.c.) 
+ ^{2{\Fll\' + \Frr\^) + \Flr\^ + \Frl\^ 

+ {B^^FIl + B^F*^R + B^FIr + B^F*^L + h.c.) + 2(^f'F£^ + F^^ + h.c. 
+ {A^FIr + AfF*^L + h.c.) - 4{A^FIl + A^F*^R + h.c.) 

-2{A^F*^L + A^Fl^ + h.c.)} 



Fll 



FlZI 



Frr 



FrZ[ 



R 



Flr 



FlZ 



R 



Frl 



FrZI 



(3.5) 



(3.6) 
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Some box form factors have been redefined to include the contributions from the Z' penguins: 

Bf ^ B^ = Bf + 2FiL, (3.7) 

^ B^ = B^ + 2F'j,^, (3.8) 

B^ ^ B^ = B^ + FiR, (3.9) 

B^ ^ B^ = B^ + F'j,L, (3.10) 



with 



^^^"IvflT' ^^"IwIT' ^^"^vflT' ^"^-^^^ 



Z' 



In our case many of the form factors are zero so the expression above will be somewhat 
simplified. 



3.1.3 /iN— >eN conversion rate 



The /i — )• e conversion process |42y43j is similar to /i — )• eee and differs only in that the lower 
part of the diagrams is coupled to quarks instead of leptons. Also, we do not have identical 
particles in the final state. We will sort the form factors as follows: 

VW^peng = -^uM [Q^Y{A{Pl + A^Pr) + m^ia^'''Q,{A^PL + A^Pr)] u^{p), 

XUg{p2)Qq7f,Vg{p3) , (3.12) 
Mzpeng = J^uM[j^'{FLPL + FRPR)]u^{p)xu^{p2h^-^^^-^Vg{ps), (3.13) 

z 

g2 + Z"^ 

A^z'peng = ^txe(pi) {i^{F'^Pl + F'rPr)\ Mp) X ^,(^2)7/. ^ ^ ^^.(ps) , (3.14) 
z 

BL 

■^box = e'^UerPLU^,{p)xUgip2h^,Vgip3). (3.15) 

We have already taken into accounto that, out of the original box form factors, only B^ is 
non-zero. The process width is [43j: 



Z 



4 



r{fi ^ e) = a5^|F(g)|2m^ \2Z{Af - A^) - {2Z + N)Bt - {Z + 2iV)5f,| . (3.16) 
The vertex form factors are as for ^ — )• eee and were given in ()3.4p . We also have defined 

Big - B^g + + (3.17) 

to include the contributions from the Z' penguins. 

The conversion rate is obtained by dividing by the muon capture rate: 

n = . (3.18) 

-L capt 

The nuclei we will consider are 22^1 and ^^jAu, whose relevant parameters are listed in 
table [ini 
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Nucleus 


Z N Zeff Fiq) Tcapt [GeV] 


iTi 


22 26 17.6 0.54 1.7 x 10"^^ 
79 118 33.5 0.16 8.6 x 10"^^ 



Table 10: Relevant input parameters for the nuclei under study. From }44j . 




II IV V VI 



Figure 1: Diagrams for /i — t- 67. 



3.2 yu, 67 



In contrast to the LHT model, there is no symmetry preventing the coupling of two standard 
particles to a heavier one. Thus, we classify the contributions to ^ — )• 67 into two types of 
diagrams (see figured]): those involving heavy X gauge bosons and those with W bosons in 
the loop. Since only dipole form factors contribute to this process, we have: 



^ M 



Fl\x + Fl\w. 



(3.19) 
(3.20) 
(3.21) 



Defining the mass ratios: 



m 



Mi' 



CO 



Ml 



(3.22) 



we find the following contribution to the dipole form factors for the X-based diagrams: 



Fl\x- 

where we have introduced 

Fx{x) 



.p7| ^ aw rn^ 
^^^'^ 16^ M2 



Mi 



2Ci - 3Cii - X ( Co + 3Ci + -Cii 



5 3x - 15x2 - 6x3 3^3 

6 12(1 - xf ^ 2(1 - x) 



Inx. 



(3.23) 

(3.24) 
(3.25) 
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The loop functions are summarized in Appendix |Al In this case 
photon. 

For the Vl^-based diagrams, we obtain 



for an on-sheh 



where 



Fm\w 



(3.26) 



Fw{x) 



Ml X 



Co + Ci 



11 



x(-7 + 5x + 8x2) x2(-2 + 3x). 



12(1 -x) 



2(1 - x)4 



(3.27) 
(3.28) 



The total dipole form factors are therefore: 



^ M 



^ _ aw rUf, ^ .He*^rif^ 



jj2 Fx{xi) + 5^ Fw{xi/uj) 



(3.29) 



It should be noted, however, that in order to single out the relevant contributions in 
W based diagrams, it is necessary to expand the loop functions with arguments Xi/u as a 
series in a; cx f^//^ to see whether the term is of leading order (w^//^) or of higher order. 
Nonetheless, once a term is accepted, we use the full expression for the numerical calculation. 
This is because we can not safely expand in orders of a; and then take a low Xi limit since the 
limits f^//^ — )• and Xj — )• do not commute. Expanding in oj spoils the low Xj behavior 
since w = is a singular point. In this limit some terms in the series are divergent for low Xj 
(for instance oc Inxj) and give abnormally large contributions in this area not corresponding 
to the physical case, and it is therefore safer to keep the full expressions for the loop functions. 
We apply these criteria to all W based contributions discussed in this paper. 



3.3 — )■ eee 

Here we have contributions from diagrams with 7, Z and Z' penguins as well as boxes. The 
relevant diagrams are displayed in figure [2l 



3.3.1 Photon penguins 

The dipole form factors are the same as in the — )• e7 case. For those terms we can set 
= as in /X — )• e7, since is small in ^ — )• eee. However, the Fl and -Fr form factors 
require loop functions up to order QP' jM"^ where M is the mass of the gauge boson in the 
loop, because these terms need to cancel the factor coming from the photon propagator. 

In and Fr we neglect terms of order m^jM'^ . This means that, Fr ~ 0. We again 
divide the contributions into two groups. Firstly, the X related diagramas work out as 

= ^ E "^r^r Gx{x.), (3.30) 

i 
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Figure 2: Relevant triangle and self-energy diagrams for /x — t- eee. Others give subleading 
contributions. 
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where 

Gx{x) 



and Ag = I 



(1 - Si) 



+ 5i + 6Coo + x[-Bi + Coo - MlCo 



A,-ln 



M' 



Ml 



5 x(12 + x-7x2) a;2(12- lOx + a;^) 
~18^ 24(1 + 12(1 ' 

7 + In 47r, divergent in four dimensions. Tlierefore: 

aw 



F? 



Air Ml 



- (^2Ci + iCii j Q2 
(3.31) 
(3.32) 



(3.33) 



The V^-based diagrams contribute with 

aw 
Ait 



Fl\w = ^ VrV^'^ Gw{x,/u) , 



(3.34) 



where 



That is: 



Gw{x) = ~ + Bi + 6Coo + x(^Bi + Coo-M^Co 



G^^ix) 



+ 



X 



F7 



l\w 



1 _ x(-2 + 7a; - llx^) 
6 72(1 - xf ' 12{l-xY 

aw 



\n.x . 



47r M^ 



(3.35) 
(3.36) 

(3.37) 



3.3.2 Z penguins 

There are three pieces in this section: two of them {F^\x and F^\w) involve only heavy neu- 
trinos in the loop, and the third (-F^|hi) contains one heavy and one light neutrino exchanging 
either a or an X boson, 



Fi = Fi\x + Fi\w + F^\u. 



(3.38) 



Again 2± if we neglect rnf/M"^. The X-based diagrams result in: 



where 



aw 1 
Att cwsw 



Ix{x) 



Ix{xi) + 5l Hx{xi) 



6x — x^ 2x + 3x^ 
2(1- x) ~ 2(1- x)2 



Inx, 



(3.39) 
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Hx{x) = l + ^^^^lnx. (3.40) 



The W boson diagrams give the following contribution: 



Pl\w = (>u 

4tt swCw 

t 

where 



TT r / \ 2 + (1 - tly)ti3 2 T t , \ 
Hw{xi/uj) 0^ Iw{xi/uj) 



Hw{x) = - + — - + — -^Inx, 3.41 

8 4(1 — x) 4(1 — xY 

2 2 

Iw{x) = -^ + —^—lnx. (3.42) 
1 — X (1 — xj^ 

Finally, diagrams where the Z couples one heavy to one light neutrino contribute with: 

= ^^—Siy^Vrv;^\Coo{M^,0;x,/oj)-Coo{Mj„0;x,)} (3.43) 

I 

aw 1 .2 



Air swcw 

where 



5l yrVt [Hz{x^M - Hz{xi)] , (3.44) 



Hzix) = . (3.45) 

2(1 — X) 



3.3.3 Z' penguins 

Here we have two contributions: 



Ft = Fi'\x+Fi'\w. (3.46) 



There is no piece analogous to -F^lhi since the Z' has an additional f^//^ from its propagator 
that makes those terms subleading. The form factors read: 

Ft\x = ^ ^=Y,Vrvrix{x.), (3.47) 

fE'Iw = ^ ^=Y,Vrvriw{x,/u:). (3.48) 

swy^-tw i 

where Ix and I\y are defined in (|3.39p and (j3.42p . 
3.3.4 Box contributions 

Only W and X particles can be involved in the loop (see figure [3]). Crossed diagrams, not 
shown in the figure, contribute a factor 2 due to Fierz identities [35]. Neglecting me/M we 
have contributions only to the form factor, divided in three terms: 

Bi = B^\x + B^\w + B^\wx , (3.49) 
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Figure 3: Box diagrams for fi — )• eee. 



where 



Bi\x 



Bi\w 



Bi\wx 



aw 


1 


1 


Stt 




Mi 


aw 


1 




Svr 






aw 


1 




Stt 







ij 



l;do{uJ, Xi, Xj) - 2do{uJ, Xi,Xj) 



and 



(3.50) 
(3.51) 

(3.52) 



(3.53) 



3.4 /i — e conversion in nuclei 



The triangle form factors are the same as in the /U — )• eee process. Only the box form factors 
need to be recalculated. These are of course embedding dependent. The diagrams include 
all combinations of quarks and gauge bosons (figure S]). We have assumed heavy quark 
degeneracy to suppress any mixing effects from this sector, which were not included in the 
Feynman rules anyway. In this approximation, only diagrams with a D quark appear in the 
anomaly-free embedding while only diagrams with a U quark are included for the universal 
embedding. Diagrams with light quarks appear in both embeddings but will be found to be 
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Figure 4: Box diagrams for /xN — )• eN. 



a subleading contribution. The form factors B^^ refer to diagramas where the q quark enters 
on the lower hne. 

In the anomaly-free embedding we obtain: 
aw 1 1 



4 + -XiXD ] do{Xi,XD) 



B 



Id 



-2XiXD do{xi,XD) 
+Su{Sd + S^)XiXD 

0, 



|2x2 



+ 



4a;2 



XiXD do{xi/u},XD/^) 



do{uJ,Xi,XD) - - do{u},Xi,XD) 



and in the universal embedding: 



B 



In 



B 



Id 



0, 

aw 1 1 
16^ 4; 



1 + -Xixu ] do{xi,xu) 



-2xiXu do{xi,xu 



2x2 



+ 



4a;2 



-XiXu doixi/Lv,xu/uj) 



1 



dQ{uj,Xi,xu) - - do{uj,Xi,xu) 



with 



XD = mD/Mx, xu = mjj/Mx- 



Here, the Z' couplings ^'j/ _R appearing in (j3.17p are also embedding dependent. 



(3.54) 
(3.55) 

(3.56) 
(3.57) 

(3.58) 

(3.59) 

(3.60) 
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4 Numerical results 



We now give the results for the different processes. We follow closely the analysis in |181ll9j 
and restrict ourselves to the case of two lepton families to illustrate the physical behavior of 
the amplitudes. This leaves us with four basic parameters: 

• The masses of the two heavy leptons, rriN^^^ parameterized in terms of an average mass 
and a relative mass splitting: 

x = ViT^, x^ = T§^ S = ^ -■ (4-1) 

• The only remaining mixing angle ^ in V^: 

(cos 9 sin 9 \ 
(4.2) 
— sin 9 cos 9 1 

• The Little Higgs breaking scale /. 

Values in the theory for these parameters which may be considered natural are / ~ 1 
TeV, X = 5 = 1 and sin 2^ = 1. We take those as reference values for our analysis as these are 
the typical values one could expect for this model. There are also two additional parameters 
involved: and the mass of the heavy quarks that appear in the box diagrams of the 
N — )• e N conversion process. We will parameterize the heavy quark masses as xjj = mfj/M'^ 
and XD = m'j^/M'^ and will take reference values xu = xd = 1- In order to comply with the 
limit 5i, < 0.05 we need to impose ftj^ > 3.5 TeV. Therefore, to allow for values of / below 1 
TeV, we choose a reference value = 5. 

For the processes — )• e7 and — e conversion, all the form factors have the following 
general form: 

^ = E VrvrFix.) = ^[F(xi) - F{X2)] . (4.3) 

1=1,2 

The branching ratio or conversion rate can then be approximated by: 

9 2 



-T:,5sm29 



(4.4) 



The dependence on sin 20 is exact while the dependence on 6 is valid for small values of 
this parameter. The dependence on / is also approximate due to the fact that some loop 
functions depend on Xi/uj. However, the leading order in every form factor is f^//^ and the 
1 //^ behavior is a good approximation. The dependence on x cannot be expressed as simply. 

For the — t- eee process, the vertex form factors are as in ()4.3p while the box form factors 
take the form: 



ij=l,2 

sin 26* 



cos^ 9 {F{m\^ , ) - F{m\^ , m\ )] + sin^ 9 {F{m\ , ) - F(m\^ , )] 

(4.5) 
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fj, eee 


fj. Au e Au 


^Ti^ eTi 


Limit 


1.2 X 10-11 


10-12 


7 X 10-13 


4.3 X 10-12 


//TeV > 


4.3 


4.1 


13.9 (16.5) 


8.8 (10.3) 


sin 29 < 


0.052 


0.055 


0.005 (0.004) 


0.013 (0.009) 


\6\< 


0.050 


0.059 


0.005 (0.003) 


0.013 (0.009) 



Table 11: Bounds on SLH parameters from present |46fl49j limits on LFV processes. For 
N — )• e N the numbers correspond to the anomaly- free (universal) embedding. 







fi — )• eee 


fiTi- 


>eTi 


Limit 


10-13 


10-14 


10" 


18 


//TeV > 


14.2 


12.9 


397 


(468) 


sin 26 < 


0.004 


0.005 


< 10-4 


(< 10-4) 


\6\< 


0.005 


0.006 


< 10-4 


(< 10-4) 



Table 12: Bounds on SLH parameters from future |3 11134] limits on LFV processes. 

This makes the angle dependence in this process slightly different from the other cases. How- 
ever, the effect is small and the general behavior remains similar to (14. 4p . 

In figure [SI we plot the branching ratios normalized by their current experimental limits as 
functions of each of the parameters /, x, 5 and 9. In each case, we vary one of the parameters 
and keep the remaining ones at the reference values. Only normalized values below the unity 
are experimentally allowed. The general behavior of the SLH model is very similar to that of 
the LHT model studied previously in |18pi9j. Tables [TT] and [T2l contain generic constraints on 
the scale /, 5 and sin 29 for standard values for the rest of the parameters in each case. As can 
be observed in the x plot, there are cancellations between different contributions to /uN — )• eN 
conversion processes, but for rather low heavy neutrino massesQ The abrupt change of the 
slope along the curves away from zero reflects the change of the sign assignment e in ()2.63p . 
which we choose to minimize the size of the conversion rate. Figure [6] shows the dependence 
on tjs. Notice that, in general, this parameter does not allow us suppress the branching ratios 
enough to get within the experimental bound. Furthermore, the dependence on is very 
light for larger values of this parameter due to the fact that it appears in the denominator of 
5^, and simply suppresses these terms for large values. Low values for makes the calculation 
unstable (the expansions of the scalar fields would no longer be good aproximations) , being 
forbidden anyway because of the limit 5i, < 0.05. 

*In this example, the cancellations occur from interference among boxes and penguins. A similar effect 
can be seen in the LHT model (see figure 4 of Ref. [IS]), where the 7 and Z penguins can equate the box 
contribution in absolute value at specific points but with opposite sign. In the same figure, one can also see that 
the penguin contributions vanish at some points, where they flip sign. Although we only show cancellations in 
N — >■ e N conversion, they are not necessarily limited to this case and it is possible that different combinations 
of parameters could generate cancellations in other processes. 
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Figure 5: Ratios of SLH predictions to current limits with = 1 and xu = xd = 1- 
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In figure [7| we show exclusion contours for current experimental limits in the {5, sin 29) 
plane. Points below the contour lines are within the measured bounds for all studied LFV 
processes although, in general, the fiAu — )• eAu process gives the most stringent limits. 
These plots show that the mass splitting and the mixing angle must be correlated in order 
to suppress the LFV effects. The correlation is similar for both embeddings, both of which 
require very small mixing angles or mass splittings to stay within experimental constraints. 



— > 67 
/i — )• eee 
/^N -> eN : 
(AF) Ti 
(U) Ti 
(AF) Au 
(U) Au 



i2 



Figure 6: Ratios of SLH predictions to current limits as functions of with natural values 
for all other parameters. The shaded region on the left is excluded by the limit on the mixing 
of SM leptons with heavy neutrino singlets. 

In figure[8]we show scatter plots for the different processes. The points are generated from 
a random scan of the parameter space which was done taking x uniform in the range [0.25, 4], 
6 uniform on the range [—4,4], uniform in a range corresponding to / G [0.5,50] TeV, 
sin 20 uniform in [0, 1] and tj^ uniform in [0.1, 10]. We then discard points that do not comply 
with the limit 6i, < 0.05, i.e., ftjs > 3.48 TeV to obtain the red points. The green points 
are generated by also excluding points that do not fulfill the current limits on the process 
which is not plotted. For instance, if we are plotting /i — t- 67 against fj, — >• eee, we exclude 
(red) points that have the /xN ^ eN process over its current experimental bound. The blue 
points are analogous to the green points but we use the future limits on the complementary 
process rather than the current limit. We show /i — t- e7, ^ — )• eee and /iN — ?■ eN on Ti (the 
future limit for this nucleus is the strongest) in the anomaly free embedding (both embeddings 
produce similar results and we choose just one for simplicity). The shaded area is allowed 
by current experimental limits on the plotted processes. No blue points are plotted in the 
first graph because the future limit on /i N — )• e N on Ti excludes practically all points in the 
scanned range. In all cases, we see that the branching ratios are clearly correlated with each 
other due to the fact that the basic behavior is the same in all cases, as given in equation 
(14. 4p . The only deviations from this come mainly from the dependence on x and, to a lesser 
degree, from tyg which is different in each case. The range chosen for this parameter roughly 
fixes the width of the scatter plot while the ranges for the other parameters (notably /) fix 
the length of the plot. For illustration, we draw the curves of constant / = 1, 10, 50 TeV and 
varying x. The thicker part of these lines correspond to the allowed x range in the scatter 
plots [0.25, 4]. The edges of the point distribution approximately align with the edges of these 
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-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 1.0 



S S S 

Figure 7: Exclusion contours in the {S, sin 2^) plane for the anomaly free embedding (top) 
and the universal embedding (bottom). Each curve corresponds to a different value of / (from 
the bottom up / = 0.5, 1, 2, 3, 4 TcV) and each plot to a different value of x (from left to 
right X = 0.5, 1, 4). Also xu = xd = 1 and = 1. 



33 



T ' 1 ' 1 ' 1 1 1 1 1 1 1 1 1 1 1 1 r 




10-18 10-16 10-14 10-12 10-10 10-8 10-6 



— > eee) 

Figure 8: Scatter plots for random scans over model parameters (see text). 
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thicker sections. The ahgnment would be exact if the scahng ()4.4p were exact and if there 
were no dependence on t^. In the cases where we plot //N — )• eN we see the effect of the 
cancellations from figure which translates into the thin lines approaching vanishing values 
of the conversion rate for this process. The point density is not representative of the relative 
frequency of each point type since they are obtained from separate runs (they are not subsets 
of a single set of points) and are forced to have a specific number of points each. Also, different 
choices of the parameter distributions can produce other point densities. The purpose of the 
plot is to illustrate the area reachable while keeping the parameters in the specified ranges. 
This means that, to some extent, experimental detection of one of these processes constrains 
the possible values for the others. 

5 Conclusions 

Although lepton flavor mixing in the SLH model had been previously addressed [5], no phe- 
nomenological studies of the basic flavor changing proceses had been presented before. To 
this end, a thorough analysis of the Lagrangian has been performed in order to obtain all the 
required field interactions in the 't Hooft-Feynman gauge, which is well suited for one-loop 
calculations. In particular, we had to identify to the desired order in v/ f the actual would- 
be-Goldstone bosons, those to be eaten by the corresponding physical gauge bosons. Then 
the Feynman rules for all needed interactions were obtained. 

We have analytically calculated the amplitudes of the basic lepton flavor changing pro- 
cesses in terms of standard one-loop functions, that have been reduced to relatively sim- 
ple expressions. These amplitudes are ultraviolet finite, as they are in the case of the 
LHT model |18 p i9 t [2 H l22j. On the other hand, the heavy neutrino Yukawa exhibits a non- 
decoupling behavior (figure [5]): both Z and Z' penguin contributions to the amplitudes grow 
with mj^/M'^, a result that is well known and was discussed in |18yi9]. 

To simplify the phenomenological analysis, we have assumed, as in previous studies, that 
only two lepton families mix. This leaves us with just four free parameters. Present limits 
on the considered rare processes translate into bounds on these parameters that allow us to 
assess the degree of naturalness of the model. From table [TT] we conclude that the value of the 
breaking scale / must be above 14 TeV, which is more stringent than the limit derived from 
EWPD, about 4 TeV [6H9]ll| Thus, near our standard reference values, flavor constraints in 
the SLH model and in the LHT model are similar. For these specific values, the strongest limit 
on the SLH / scale is of almost 14 TeV for the anomaly free embedding while the analogous 
limit in the LHT model is of about 10 TeV [l9j. The universal embedding is similar with a 
limit of about 16 TeV. However, we must keep in mind that these limits are not rigid and 
depend on the reference point. Moreover, the point was different for LHT than for SLH {x 
was taken at 4 in the LHT in order to avoid an unnatural cancellation). Additionally, for 
equal /, the gauge bosons in the LHT model are heavier than those in the SLH model by 
about a factor ^/2. Keeping these differences in mind we find that the results for either model 
are similar and of the same order of magnitude. 

t This bound also applies to variations of the original SLH model '25| implementing a lighter new T quark 
and a richer Higgs phenomenology consistent with present data |50ii51) . Such an extension does not alter the 
flavor structure discussed in this work. 
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As can be observed in the scatter plots in figure [HI present LFV experimental limits seem 
to allow a relatively large region of the parameter space (green and blue points in the shaded 
regions), as is the case for the LHT model |17y22j. This does not contradict our main claim 
about the necessary fine tuning. Indeed, in order to be in the experimentally allowed area the 
effective misalignment between the light and heavy lepton flavors 5 sin 20, depending not only 
on the rotation angle but also on the heavy mass difference, must be very small and at the per 
cent level (or otherwise / of the order of several TeV and relatively large). On the other hand, 
there is a large correlation in both models between fi ^ ej and fj, — t- eee predictions. This is so 
because both processes scale to a very good approximation as \jjS sm26\'^ (see equation I4.4p 
and the dependence on the other parameters is mild within their natural range. The situation 
is similar for fj, — e conversion because in the natural region there are no large cancellations 
among different contributions. However, in general there are narrow regions in parameter 
space where these processes can vanish as illustrated in figure [5] and in figure [8] by the black 
curves going to zero values of the conversion rate for this process. Obviously, there are many 
other parameters in a formulation with three families but the behavior will be analogous using 
a convenient parametrization for the misalignment flB]. The former constraints significantly 
restrict model building, giving rise to a little flavor hierarchy problem also in LH scenarios. 
The fact that these models are so sensitive to present LFV processes, however, also implies 
that they could explain possible future observations of LFV by MEG [3T1[32] at PSI or by 
PRISM/PRIME [MllM] at J-PARC. 

Finally, although for simplicity we have neglected any flavor violation in the quark sector, 
this aspect deserves special attention. In particular, it may explain a possible discrepancy 
between the SM prediction and the measured value of quark rare processes, like for instance 
the muonic B decays — )• fip, and 5° — )• K* fiji [52] (see [15] for a review and further 
references), to be precisely measured at LHCb [MllSl]. Present upper bounds on the former 
by the CDF [55] and DO [56] Collaborations are still almost one order of magnitude above 
the SM prediction, whereas the asymmetry measurements of the latter by BABAR [57j , 
BELLE [58] and CDF |59ll60j may hide a hint of new physics. At any rate, the quark sector 
of the SLH model allows for a quite rich phenomenology because in contrast with the LHT 
model, there are FCNC already at tree level, only suppressed by small mixing angles, and 
eventually of the same size as the one-loop effects resulting from the exchange of new particles 
only. As we have discussed, one can consider two different completions of the strong sector, 
one universal with an extra vector- like quark of charge 2/3 per family and an anomaly free 
completion with an extra vector-like quark of charge 1 /3 for each of the first two families and 
an extra vector-like quark of charge 2/3 for the third one. In the universal case there are 
no tree-level flavor changing Z couplings in the down sector but there are in the up sector, 
whereas there are tree-level FCNC for both types of quarks in the anomaly free case |6HI62]. 
In either case there are one-loop FCNC also contributing to the quark transitions. Studying 
the corresponding flavor constraints would be also interesting in order to discriminate this 
model from other SM extensions [35il36j . Let us note too that although the LFV bounds 
may be a priori more stringent than those from quark flavor violation, the SLH contributions 
to charged lepton transitions are one-loop suppressed and the corresponding limits on some 
model parameters may be eventually comparable to those derived from hadronic tree-level 
processes. 
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A Loop functions 



The most general 3-point functions can be written as C{p1,Q'^,P2',Mf,M2,M^). In our 
case Pi = P2 = and only the following general types are relevant for our analysis: C = 
C(0,Q^0,M2,M|,M|), C = C{0,Q^,0,Ml,Mf,Mf) and C = C(0, 0, Mf, M|, 0) (the 
last one, symmetric under the exchange Mi -H- M2). We define the mass ratio x = M|/M^ 
and reparameterize the functions in terms of M^, and x. The resulting functions are the 
same as those in [I8l[l9]. We recall their expressions: 



Co(M2,Q2.^) = _L 



+ 



1 — X -|- In X 

(1 -X)2 

Q2 -2-3x + 6x^ 



6x Inx 



Co(M2,Q2.^) = _ 



+ 



Ml 
1 



I2x{l-xf 
— l + x — x\nx 

(1 -X)2 

4 + 6x - 3x2 _ _^ g^2 
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(A.l) 

(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 



(A. 



The 2-point functions are written in a similar way. Their general form is i?(p2;M2,M|), 
but we only need functions B = B{0; Ml,M^) and 5 = 5(0; Mf, Ml). Only Bi IS necessary. 

, M?\ 3-4x-Fx2-K2x2lnx , , 

^^7^J+ 4(13^^^ • (^-^^ 



5i(M2;x) 



-1 



A, 
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The 4-point functions are given by D{pl,pl,pl,pl, (^1+^2)^,(^2+^3)^; Mq, , M^^Ml) 
and we need only functions with zero external momenta Z?(0, 0, 0, 0, 0, 0; Mq , M^, M|, M|). 
These functions are then symmetric under any exchange of masses. We define do = MqDq and 
do = AMqDqq so that these functions now depend only on 3 mass fractions: s/^^o which 
we denote x, y and z respectively. The general expressions required are as follows [TS1[IU]: 



do{x,y,z) 



do{x,y,z) 



xlnx 



ylny 
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The case in which z 
do{x,y) = 
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1 is also necessary: 
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